Time-dependent energy absorption changes during ultrafast lattice deformation 



O 

o 

(N 



H. 0. Jeschke, M. E. Garcia and K. H. Bennemann 
Institut fur Theoretische Physik der Freien Universitat Berlin, Arnimallee 14, 14195 Berlin, Germany, 

(February 6, 2008) 

The ultrafast time-dependence of the energy absorption of covalent solids upon excitation with 
femtosecond laser pulses is theoretically analyzed. We use a microscopic theory to describe laser 
induced structural changes and their influence on the electronic properties. We show that from the 
time evolution of the energy absorbed by the system important information on the electronic and 
atomic structure during ultrafast phase transitions can be gained. Our results reflect how structural 
changes affect the capability of the system to absorb external energy. 
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I. INTRODUCTION. 



For systems in thermodynamical equilibrium the oc- 
currence of phase transitions can be inferred from a 
change in the behavior of the thermodynamic potentials 
(or its derivatives) as a function of an intensive vari- 
able. For instance, divergences or abrupt changes in 
the heat capacity C = dU/dT are indicators for most 
first-order phase transitions. In systems subject to ex- 
treme nonequilibrium conditions like excitation with an 
intense femtosecond laser, neither thermodynamic func- 
tions nor intensive variables are well defined. However, a 
large variety of nonequilibrium phase transitions can be 
induced by femtosecond laser pulses, as has been demon- 
strated both experimentally and theoretically [|l]-§| . Such 
ultrafast phase transitions are detected indirectly by 
time-dependent changes of the optical properties (j-g] 
or directly by post mortem analysis M and recently 
through more direct measurements using ultrashort X- 
ray pulses ||^] . In a recent experimental study, Musella 
et al. reported results on the excitation of graphite with 
very long pulses (in the millisecond range). The thermo- 
gram recorded during the heating of the sample shows 
a sharp increase of the temperature slope at the melt- 
ing point, which corresponds to changes in the energy 
absorption U a bs(t) and heat conductivity upon melting. 
Note that in these experiments, thermal equilibrium is 
reached during the excitation. Note that in the case 
of femtosecond excitation E a ^ s (t) is a physical quantity, 
whereas C and T are no longer well defined. The ques- 
tion which immediately arises from this analysis is: Does 
E^bs (t) show similar changes as E/ a b s in the case of ther- 
mal equilibrium when the system is illuminated by a fem- 
tosecond laser pulse? 

In this paper we show that this is the case, and that 
-^abs exhibits dramatic changes as a function of laser in- 
tensity I when the system undergoes an ultrafast phase 
transition. This leads to a natural generalization of "heat 
capacity" for laser induced nonequilibrium processes as 
C = dE ahs /dI. 

In Sec. p we briefly review our method of calculat- 
ing the structural response of covalent materials during 
and after intense laser excitation. In Sec. Ill we present 



results for diamond and graphite, concentrating on the 
efficiency with which energy is absorbed from the laser 
pulse. We show that the structural changes during the 
action of the pulse are particularly important for the un- 
derstanding of the nonlinear dependence of the absorbed 
energy on the intensity of the pulse. As these structural 
changes occur on a time scale of approximately 100 fs the 
pulse duration is responsible for very different absorption 
behavior in the r = 20 fs to r = 500 fs pulse duration 
range. Finally, in Sec. IV we summarize our results. 



II. THEORY. 

The calculation of the energy absorption of covalent 
solids is based upon a theory for the analysis of laser in- 
duced ultrafast processes in solids. The theory employs 
a molecular dynamics method on the basis of an elec- 
tronic tight-binding Hamiltonian. This real-space calcu- 
lation takes into account all atomic degrees of freedom. 
Special attention is paid to the strong nonequilibrium 
created in the electronic system by the ultrashort laser 
pulse. A method of calculating nonequilibrium occupa- 
tion numbers for the energy levels of the system leads 
to a molecular dynamics calculation on time-dependent 
potential energy surfaces. This approach provides a theo- 
retical framework for the treatment of strong nonequilib- 
rium situations in materials where atomic and electronic 
degrees of freedom play an equally important role. 

In this theory trajectories for laser induced structural 
changes are determined by molecular dynamics simula- 
tions for which we calculate forces from the Hamiltonian 
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The second term contains a repulsive potential 

that takes care of the repulsion between the ionic cores, 

and the first term is a tight-binding Hamiltonian 
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Here, the first term is the on-site contribution. In the 
second term, the & are the hopping integrals, and cf 



and c -^ are creation and annihilation operators for an 
electron at site i or j in orbitals r) or i9, respectively. 

From the Hamiltonian, Eq. p) , we calculate forces us- 
ing the HcHman-Feynman theorem: 

ffc({ry(t)},t) = -^n(e m ,t)(m|V fe iJ T B({ry(t)})|m) 
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diffusion processes is assumed: 
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Here, m) stands for the eigenvector of the Hamiltonian 
.Htb that corresponds to the eigenvalue e m . The special 
feature of the theory which makes it applicable to opti- 
cally excited materials is contained in the time dependent 
occupation numbers n(e m ,t) for the energy levels e m of 
the system. While in thermal equilibrium these occupa- 
tion numbers are calculated from a Fermi-Dirac distribu- 
tion function n°(e m ) — 2/(1 + exp{(e m — /«)/£;b7c}) at 
a given electronic temperature T e , electronic nonequilib- 
rium is accounted for by solving equations of motion for 
the occupation of electronic states: 
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dt 
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Thus, the electronic distribution is at each time step 
folded with the pulse intensity function g(uj,t). This 
means that at each time step, the occupation of an energy 
level e m changes in proportion to the occupation differ- 
ence with respect to levels at e m — fuv and at e rn + hw. 
In Eq. (0), constant optical matrix elements are assumed. 
The second term of Eq. fl) describes the electron-electron 
collisions that lead to an equilibration of the electronic 
system with a rate equation of the Boltzmann type for 
the distribution n(e m , t). Hence, with a time constant n, 
the distribution n(e m , t) approaches a Fermi-Dirac distri- 
bution n°(e m ). 

The electron temperature T c that results from the elec- 
tron thermalization will not remain constant over time, 
but will decrease due to electron-phonon coupling and 
to diffusion of hot electrons out of the laser excited re- 
gion of the solid into colder areas. In the regime of an 
excitation of typically 10 % of the valence electrons into 
the conduction band for which this theory is intended 
no precise knowledge about thermal conductivity of the 
electrons and electron-phonon coupling constants is avail- 
able. Thus, a relaxation time ti which characterizes the 
time scale of the total decrease of the electron tempera- 
ture by both the electron-phonon and by the hot electron 
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This simple relaxation time approach has the advantage 
that the time scale ti on which the electron-lattice equi- 
libration takes place can be taken from experiment. 

The forces given by Eq. (g) can now be used to solve 
the equations of motion for the atoms numerically. In 
the case of a bulk system, a constant pressure molec- 
ular dynamics (MD) scheme is used. This is based on 
a Lagrangian which contains the shape and size of the 
molecular dynamics supercell as additional degrees of 
freedom JL0|: 
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The first term is the kinetic energy of the atoms, with 
the coordinates of the atoms r; = hsi written in terms 
of the relative coordinates Sj and the 3 x 3-matrix h that 
contains the vectors spanning the MD supercell; Sj are 
the relative velocity vectors, and T denotes transposi- 
tion. &({rij},t) is the potential which is calculated from 
a tight-binding formalism as explained above. The terms 
in the Lagrangian of Eq. (ra) which are responsible for 
the simulation of constant pressure are an additional ki- 
netic energy term i^ ce ii for which the simplest form is 
K C ci\ = ^z 1 Ti(h T h) fll0| , and U cc n is an additional po- 
tential term which describes the effect of an isotropic ex- 
ternal pressure U cc u — p cx t^; ^ = det(/i) is the volume of 
the MD supercell. The equations of motion are derived 
from Eq. (g) by the Euler-Lagrange formalism, and they 
are integrated numerically with the velocity form of the 
Vcrlct algorithm |I|@. 



III. RESULTS. 

The dependence of the absorption of energy on the 
laser pulse intensity and duration shows an interesting 
behaviour. The following results have been calculated 
in the bulk and for constant pressure. Diamond and 
graphite samples of N — 64 atoms per MD supercell were 
thermalized to a temperature of T = 300 K by simulated 
annealing before absorbing laser pulses of a wide range 
of intensities and durations. 

Results for absorption of carbon in its diamond struc- 
ture are shown in Fig. nl. Every square (□) corresponds 
to a calculated trajectory. Note, the lines connecting 
the dots are drawn to guide the eye. The graph rep- 
resents the actual absorbed energy in units of electron 
volts per atom as a function of the laser intensity for 
different pulse durations. This quantity is denoted here 
as the "offered" energy. It is the energy that would have 
been absorbed by the system, if no structural changes had 



taken place. Thus, the "offered" energy is determined by 
performing the calculation of energy absorption for fixed 
atomic coordinates. Note, structural distortions change 
the electronic structure of the material and consequently 
the absorption characteristics. 




FIG. 1. Absorbed energy per atom as a function of pulse 
intensity ("offered" energy) for diamond. The pulse duration 
ranges from r = 20 fs to r = 500 fs. The shape of the ab- 
sorption curves and their ordering with respect to the pulse 
duration can be understood from the time-dependent changes 
of the DOS. Each of the results (□) corresponds to a calcu- 
lated trajectory. 



A general trend shared by all but the r = 390 fs to 
t = 500 fs curves in Fig. [j] is a moderate, linear increase 
of absorption with intensity up to an absorbed energy 
of roughly E a — 1 eV/atom. Then a sharp increase of 
absorption until approximately E = 12 eV/atom oc- 
curs. Then we find again a moderate linear increase of 
absorption with laser intensity. The results show that the 
absorption increases more sharply with increasing pulse 
duration. With increasing duration of the laser pulse, the 
laser intensity for which the sharp increase in absorption 
occurs decreases significantly. The curves for long pulse 
durations are even characterized by the absence of an 
initial weak increase of absorption with intensity. 



In Fig. we show the absorbed energy of diamond as a 
function of time. The laser pulse duration was r = 350 fs. 
This figure serves as an example to clarify how the form of 
the absorbed energy curves of Fig. p] comes about. Fig. 
is directly related to the r = 350 fs curve of Fig. [fl: Each 
final absorbed energy at t = 700 fs in Fig. corresponds 
to a result (□) in the r = 350 fs curve of Fig. |l| 
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FIG. 2. Absorbed energy per atom in diamond as a func- 
tion of time for a pulse duration of r = 350 fs. Absorption of 
trajectories with different final states are shown. The dashed 
lines correspond to the absorption that is expected if no struc- 
tural changes occur during the action of the pulse. 



Fig. @ shows the absorbed energy for four different in- 
tensities. The final state of the material is indicated. The 
lowest pulse intensity does not lead to structure changes, 
due to the two medium intense pulses, the diamond struc- 
ture is damaged, and the strongest pulse leads to ablation 
of the material. In Fig. [| the absorption that would have 
been expected if no structural changes had occurred dur- 
ing the action of the laser pulse is indicated by dashed 
lines. We observe that the deviation of the energy actu- 
ally absorbed by the material from the expected value in- 
creases with increasing pulse intensity. In the case of the 
trajectory where diamond remains unchanged, the actual 
absorption is very close to the expected value. However, 
in trajectories where structural damage or even ablation 
is observed the actual absorption sharply deviates from 
the expected absorption. 

The cause for this sharp deviation will be explained 
with the help of the densities of states shown in Fig. [|. 
They correspond to the curve of Fig. |] with the second- 
highest absorbed energy. In Fig. g (a) the DOS before 
the laser pulse at t = —700 fs is shown; the time is mea- 
sured with respect to the maximum of the laser pulse 
of t = 350 fs duration. The DOS corresponds to dia- 
mond at T = 300 K. The DOS in Fig. | (b) is taken at 
t = 120 fs which corresponds to the time of the sharp 
bend in Fig. 0. In the DOS, the emergence of states in 
the Eg = 5.5 eV gap of diamond can be observed. This 
explains the strong increase in absorption. Fig. (c) 
shows the DOS at t = 240 fs. At this time the inten- 
sity of the laser pulse has already decreased significantly 
from its maximum value, but in the DOS, states have 
completely filled the diamond gap, and the absorption of 
the material is far stronger than before the pulse. 
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FIG. 3. Density of states for diamond at different times 
for a pulse duration of r — 350 fs. This figure corresponds 
to the trajectory with the second highest absorbed energy of 
Fig. H. For a discussion of the features of the DOS, see the 
text. 



An explanation of the overall behaviour of the ab- 
sorption in the material must be given by the time- 
dependent changes in the density of states. Note, this 
is the main material property taken into account in the 
treatment of the absorption. This time dependency of 
the DOS is caused by the structural changes induced by 
the laser pulse. Inspection of the trajectories shows that 
the sharp bend of the graphs close to an absorbed energy 
of Eq = 1 eV corresponds to the damage threshold of the 
material. The damage (disorder) causes the E g = 5.5 eV 
gap present in diamond to be filled with states until the 
gap-less DOS is reached. For pulses that are long enough 
to permit the structural changes to occur before the in- 
tensity has fallen off completely, the electronic states that 
fill the gap of diamond will enormously increase the ab- 
sorption of the material. Thus, the r = 20 fs graph shows 
only a slight increase of the absorption above the damage 
threshold, because the pulse has already run out before 
essential changes in the structure have taken place. On 
the other hand, the r = 390 fs to r = 500 fs graphs 
show an absorption growing sharply with pulse intensity, 
because all structural changes induced by the laser can 



take place while the pulse intensity is still high. Thus, 
the changed DOS has a strong influence on the energy 
that can finally be absorbed from the pulse. 

In Fig. [|, the absorption behaviour for a graphite sam- 
ple is shown. The data are presented in the same way as 
in Fig. [lj. For graphite, the absorption for a given pulse 
duration depends on the pulse intensity in a fashion simi- 
lar to that found in diamond. Again, a strong increase of 
absorption takes place at the damage threshold, which in 
graphite is found to be at roughly 3 eV. At this threshold, 
melting of graphite and the destruction of the graphite 
planes takes place. At lower intensities, absorption shows 
a moderate, linear increase with laser intensity. However, 
an important difference can be observed in the ordering 
of the curves. In the case of graphite, the slope of the ab- 
sorption curves is lower for higher pulse durations. Inter- 
preting this with the help of the density of states, we find 
that structural relaxation taking place during the pulse 
decreases the density of states close to the Fermi level 
and thus prevents absorption of energy from the laser 
pulse. A physical picture for this effect can be derived 
from the bonding character: Vibrational excitation of the 
graphite planes causes atoms to be displaced perpendic- 
ular to the plane, and thus the 3-dimensional distortion 
of the 2-dimensional graphite structure causes the pure 
sp 2 bonding of graphite to be modified with an addition 
of sp 3 bonding. 
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FIG. 4. Absorbed energy per atom as a function of pulse 
intensity for graphite. This is shown for a range of pulse du- 
rations from r = 20 fs to r = 500 fs. Note, the ordering of 
the curves with respect to pulse duration shows a behaviour 
different from diamond, which reflects a different time devel- 
opment of the DOS during the pulse. Each of the results (□) 
corresponds to a calculated trajectory. 

In Fig. p^ we show the absorbed energy of graphite as 
a function of time. The pulse duration was r = 350 fs. 
Three low laser intensities produced no change in the 
graphite structure, while a fourth pulse leads to abla- 



tion of the material. For comparison, the absorbed en- 
ergy that would be expected if no structural change took 
place during the action of the pulse is shown with dashed 
lines. As in the case of diamond (compare Fig. g) with 
increasing pulse intensity the deviation of the actual ab- 
sorbed energy from the expected value becomes more 
pronounced. 
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FIG. 5. Absorbed energy per atom in graphite as a func- 
tion of time for a pulse duration of r = 350 fs. Absorption of 
trajectories with different final states are shown. The dashed 
lines correspond to the absorption that is expected if no struc- 
tural changes occur during the action of the pulse. 

In order to analyze the cause of the increased absorp- 
tion of energy we show in Fig. @ selected densities of 
states during the action of the laser pulse. These densi- 
ties of states correspond to the trajectory with the high- 
est absorbed energy in Fig. ||. At t = —700 fs, i. e. 
700 fs before the laser pulse maximum, the DOS is that 
of graphite at T = 300 K. At t = 100 fs the overall width 
of the DOS has decreased while the DOS close to the 
Fermi level has increased. At t = 200 fs, a very high 
DOS around the Fermi level has developed, indicating 
a metallic liquid phase of the material. This explains 
the strong increase of absorption during the action of the 
laser pulse. 

In Fig. |7j we show a consequence of the behaviour of 
the absorbed energy for the spatial absorption profile of 
diamond. Typically, a laser pulse has a Gaussian profile 
in space. The different absorbed energies that result as 
a consequence of structural changes in the material dur- 
ing the action of the pulse leads to the spatial absorption 
profile shown in Fig. |t]. It is remarkable that the profile 
of the absorbed energy is much steeper than the profile 
of the pulse. This may be of importance in the interpre- 
tation of ablation experiments. 

In the introduction we have suggested characterizing 
the energy absorption behavior of a system in strong 
nonequilibrium by a generalized "heat capacity" defined 



as C — dE^bs/dl. We illustrate this concept in Fig. g 
in the case of diamond. The derivative of the absorbed 
energy with respect to "offered" energy is plotted as a 
function of "offered" energy. As explained above we 
define the "offered" energy as the energy that would 
be absorbed by the system if no structural changes oc- 
curred during the action of the pulse. Thus, the "offered" 
energy is proportional to the laser intensity / and the 
quantity dE^/dEos shown in Fig. || is proportional to 
C = dE^bs/dl. The behavior of this quantity shows a 
peak which is most pronounced for long pulse durations 
(see Fig. ||). The peak position corresponds to the laser 
intensity which induces changes in the structure and in 
the DOS of the material that are most suitable to enhance 
the absorption of energy. To the left of the peak, the laser 
is not strong enough to induce structural changes in the 
material while the pulse is still active. To the right of the 
peak, the laser intensity is so high that evaporation of the 
material sets in during the action of the pulse which leads 
to decreased absorption of energy. 
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FIG. 6. Density of states for graphite at different times 
for a pulse duration of T = 350 fs. This figure corresponds 
to the trajectory with the highest absorbed energy of Fig. H. 
For a discussion of the features of the DOS, see the text. 
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FIG. 7. Spatial absorption profile of diamond for a pulse 
duration of r = 350 fs. The spatial profile of the "offered" 
laser energy is typically given by a Gaussian (dashed line). 
Considering the relation between "offered" energy and ab- 
sorbed energy in diamond for laser pulses of t = 350 fs du- 
ration (see Fig. hi, we find an absorption profile that differs 
significantly from a Gaussian (solid line). 
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states of the materials at different times during the ac- 
tion of the laser pulse shows that an enhanced absorption 
of energy is due to ultrafast changes in the DOS which 
in turn are caused by the structural relaxation of the ex- 
cited material setting in during the action of the exciting 
pulse. As these structural changes occur on a time scale 
of approximately 100 fs the pulse duration plays an im- 
portant role in determining the influence of the structural 
relaxation on the absorption behavior. 

In summary, the results show the dependence of ab- 
sorption on the light pulse form and the occurring struc- 
tural changes of these materials. These results may also 
depend on the frequency and polarisation of the laser 
pulse. Nevertheless, the fact that nonthermal processes 
may occur during the action of the laser pulse and thus 
change the absorption of the material is of general valid- 
ity. 

If we compare the theoretical result for the dependence 
of absorbed energy on pulse intensity for graphite (see 
Fig. 0) with the experimental result of Ronchi et al. |13J ] 
we find a surprising similarity. Although the pulse du- 
ration r = 20 — 30 ms in the experiment was far higher 
than that considered in the calculation, the increase of 
absorption at the melting point should have a similar ex- 
planation as the increase of absorption at the damage 
threshold in the calculation. It is remarkable that the 
effect reported by Ronchi et al. is also present for ultra- 
short laser pulses. 
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FIG. 8. Derivative of absorbed energy with respect to 
"offered" energy in the case of diamond for different pulse 
durations. The "offered" energy is proportional to the laser 
pulse intensity. The functions in the graph are calculated as 
derivatives of the functions connecting the calculated absorp- 
tion values of Fig. fil Note that the functions corresponding 
to all pulse durations between r = 20 fs and r = 380 fs have 
been multiplied by a factor 30 for clarity. 



IV. SUMMARY. 

We have calculated the absorption of energy from in- 
tense laser pulses in diamond and graphite. An inter- 
esting nonlinear dependence of the absorbed energy on 
the laser intensity is observed. Analysis of the density of 
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